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Abstract 

This paper presents a formulation of the longitudinal glideslope tracking of a transport-class 
aircraft in severe wind shear and turbulence for application to robust control system design. 
Mathematical wind shear models are incorporated into the vehicle mathematical model, and wind 
turbulence is modeled as an input disturbance signal. For this problem formulation, the horizontal 
and vertical wind shear gradients are treated as real uncertain parameters that vary over an entire 
wind shear profile. The primary objective of the paper is to examine the formulation of this 
problem into an appropriate design format for use in p-synthesis control system design. 

1. Introduction 

Wind shear is the local variation, at a particular time, of wind velocity with distance, and 
can be characterized by large rapidly changing variations in wind velocity over short distances. 
Wind shear can occur from a variety of sources, including convective outflows (such as 
downbursts associated with thunderstorms and other convective clouds), various types of wind 
fronts (such as gust fronts, sea-breeze fronts, and air-mass fronts), the jet stream, and ten-ain- 
induced wind variability. The magnitude and relative direction of a wind shear are quantified by 
dividing the velocity difference at two points by the distance between them, and this quantity is 
referred to as the wind shear gradient. When wind shear occurs at low altitudes in the vicinity of 
airport runways, it can pose a devastating hazard to aircraft in both the takeoff and landing phases 
of flight, and various documents are available which report on this flight hazard. One such report 
is the Federal Aviation Administration report by Shrager (1977), ref. [9], which reported on 
aircraft accidents or incidents related to low-altitude wind shear occurring from 1964 to 1975. A 
later report under the auspices of the National Research Council entitled "Low-Altitude Wind Shear 
and Its Hazard to Aviation" (see Committee ... (1983), ref. [5]) presents details on wind shear 
arising from all of the above wind shear sources, as well as the hazard each poses to aviation. A 
conclusion of this report is that the greatest wind shear hazard to aircraft is posed by "downdrafts 
and outflows produced by convective storms". Thus, a particularly hazardous form of wind shear 
is the downburst mentioned above, and the microburst, which is a downburst that is small in size 
(less than 2.5 miles in outflow size) and short in duration (peak winds lasting only 2-5 minutes). 
The report also states, however, that "serious aircraft accidents have also been caused by terrain- 
induced and frontal wind shears". Thus downbursts and microbursts are not the only forms of 
wind shear that are hazardous to aircraft. There are also numerous papers which report on the 
problem of detecting wind shear, such as a paper by Targ et al. (1991), ref. [11], which presents 
recent results on a remote wind shear detection sensor under development at the NASA Langley 
Research Center. This last report sites 26 transport accidents and 3 transport incidents occurring 
from 1964 to 1985 (for a combined total of 626 fatalities and 236 injuries) that have all been 
attributed to wind shear. One particularly devastating aircraft landing accident that was attributed to 
microburst wind shear occurred in 1975 at the John F. Kennedy International Airport and resulted 
in 1 12 fatalities. Another devastating encounter with microburst wind shear was a takeoff accident 
that occurred in 1982 at the New Orleans International Airport and involved 153 fatalities. The 
most recent catastrophic wind shear related transport accident occurred in 1985 at the Dallas/Fort 



Worth Airport. This accident also occurred during landing and resulted in an excess of 100 
fatalities. 

This paper presents a problem formulation for landing a transport aircraft through severe 
wind shear and turbulence using advanced robust control techniques. Although the paper by 
Belcastro & Ostroff (1985), ref. [2], presents an advanced control method for this problem, no 
specific robustness to wind shear was incorporated into the design. In addition, several ad hoc 
methods were required in the implementation of the control law during simulation in order to 
achieve acceptable results. It is postulated that the use of advanced robust control techniques 
which allow the treatment of wind shear as uncertainties could provide an improved method for 
designing control systems which are inherently robust to wind shear. However, use of these 
methods is not yet widespread, and formulating the problem for use with these methods is not 
commonly understood. This paper therefore focuses exclusively on the details of formulating the 
problem for use with advanced robust control methods. In particular, the wind shear gradients in 
the mathematical model of wind shear are treated as real uncertain parameters. The wind profile 
that will be used herein is based on a reconstruction of the wind profile that occurred during the 
Kennedy accident mentioned above. The problem will be formulated such that ^.-synthesis can be 
used to design a control law which is robust to an entire wind shear profile. A review of p- 
synthesis (and H°° theory, as well) is given in Bibel and Stalford (1991), ref. [4]. 

2. Problem Formulation 

The general uncertain system model used in this problem formulation is shown below in 
Figure 1: 




Figure 1 . Block Diagram of General Uncertain System 

where G(s) is the generalized plant, and K(s) is the controller. The vector w contains all 
exogenous inputs to the system (including commands, disturbances, and noise), u is the control 
input vector, z contains the controlled output variables (e.g., tracking errors, control position and 
rate, etc.), and y represents the measurement vector. The matrix A represents a block diagonal 
matrix of system uncertainties, and the vectors q and p are the inputs and outputs, respectively, to 
the uncertainty matrix. The nominal G(s) system is represented in state-space form by the matrices 
A, Bj, B 2 , Ci, C 2 , Du, D 12 , D 21 , and D 22 , and the additional B, C, and D matrices in G(s), 
namely B xp , C qx , D™* D qw , D qu , D zp , and D yp , represent the interconnection of the uncertainties 
into the system (see Belcastro, Chang, & Fischl (1991), ref. 4). 

The generalized plant, G(s), contains the open-loop plant, P(s), as well as any weighting 
functions used in the design. The plant, P(s), for this study consists of the aircraft and wind shear 
dynamic equations. For this problem, G(s) will be formulated as shown in Figure 2: 
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G(s) 

Figure 2. Block Diagram of General Model Subsystems 


where W w (s) represents the exogenous input weighting matrix, Asc is a constant scaling matrix 
associated with the uncertainty model, W z (s) is the performance weighting matrix, and Ga(s) 
represents the control actuator model. The overall objective of this study is to track the glideslope, 
despite severe wind shear and turbulence conditions, during the descent to landing of the aircraft. 
The subsystem models (depicted in Figure 2) as well as the associated generalized system model 
(depicted in Figure 1) required to achieve this objective within the framework of this robust control 
problem formulation are therefore developed in the following subsections. 

2.1 Plant Model Formulation 

The glideslope tracking problem was formulated for the Transport Systems Research 
Vehicle (TSRV) B-737 aircraft located at NASA Langley Research Center (see Staff ... (1980), 
ref. [10]). The linearized mathematical model of the aircraft was based on the following trim 
conditions: 


glideslope: 

weight: 

center of gravity (c.g.): 

flaps: 

airspeed: 


3 degrees 
80,000 lbs 

25% of the MAC (mean aerodynamic chord) 
40 degrees 

125 knots (211.5937 ft/sec) 


Plant Model 

The state equations for the basic plant model, P(s), represent a linearized perturbation 
model from the aircraft trim point (defined above) and are given by the following: 


where: 


x — Ax + Bu *4* B w d w 

(1) 

y - Cx + Du + D w d w + 

(2) 

x = [uwq0u w w w h]' 

(3) 

u = [ Sb Sr ] T 

(4) 

d w = l u w , w W( | T 

(5) 

y = [ V a q 0 h h ] T 

(6) 

il = [ Tlv. ilq % ilh tlh l T 

(7) 
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The first four states in the x-vector (i.e., u, w, q, and 0) represent the four longitudinal aircraft 
states in the body axes (i.e., longitudinal velocity in the x and z directions of the body axes, pitch 
rate, and pitch attitude, respectively), u w and w w represent the longitudinal and vertical wind 
velocity components expressed in the positive longitudinal and vertical directions (respectively) of 
the body axes, and h represents the aircraft altitude along the glideslope. The wind model is 
comprised of wind velocity components due to both wind shear and turbulence. The wind shear 
model describes a wind velocity vector (in the inertial Earth-fixed frame of reference) whose 
horizontal and vertical wind velocity components vary linearly with respect to decreasing altitude 
along die glideslope. Horizontal and vertical wind shear gradients, U z and W z (respectively), were 
used to represent the linear wind velocity variation relative to the aircraft's descent along the 
glideslope. The wind shear model was then rotated to the aircraft body axes and expressed as a 
function of the wind shear gradients U z and W z , as well as the associated aircraft states. A full 
derivation of this model, including the wind shear state equations, is provided in Belcastro & 
Ostroff (1985), ref. [2]. The wind turbulence model was treated as an input disturbance vector, 
d w> whose components consist of the horizontal and vertical wind gust components, u Wg and w w „ 
respectively, in the body axes. The Dryden turbulence model was used to provide a description of 
the turbulence spectrum, and this model is discussed in further detail in Section 2.3.3. The control 
input vector, u, to the plant consists of elevator deflection (8 e ) and engine thrust (87). The 
measurement vector, y, consists of aircraft airspeed, V a , pitch rate, q, pitch attitude, 0, aircraft 

altitude, h, and altitude rate, h, respectively. The vector q represents measurement noise. The 
matrix elements as well as the trim values and units associated with the vector elements for the 
above state model are given in the Appendix. 

The A and B w matrices for the above model are of particular interest in this study, because 
they contain the wind shear gradients, U z and W z , which will be treated in this problem 
formulation as uncertain plant parameters. These matrices are therefore repeated below for 
convenience: 


A = 


Aaa A aw 

A A 
''wa n ww 

A h 0 


0 

0 

0 


(7x7) 


-an 


-ai2 


B w — 


■U71 


-a72 


( 8 ) 


(9) 


The sub-blocks of A are all constant matrices except A wa (2x4), whose elements are linear (affine) 
functions of the horizontal and vertical wind shear gradients U z and W z , respectively, as shown 

below: 


a w n = .028944 U z + .00083845 W z , a w2 ] =-.00083845 U z + .028944 W z 
a w i2 = .99916 U z + .028944 W z , a w22 = -.028944 U z + .99916 W z 

awl3 = 0 , a w2 3 = 0 (10) 

a w ] 4 =-210.89 U z - 17.188 W z , a w 24 = 17.188 U z - 210.89 W z 

The A and Bw matrices are therefore linear (affine) matrix functions of the horizontal and vertical 
wind shear gradients U z and W z , respectively. The control system is usually designed for the 
nominal case of no wind shear (i.e., U z =0 and W z =0) so that the design is not "tuned" to a 
particular wind shear gradient [see Belcastro & Ostroff (1985), ref. 2]. Subsequent to designing 
the controller, these gradients are typically set to nonzero values for a closed-loop analysis under 
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various wind shear gradient combinations. For this study, however, the horizontal and vertical 
wind shear gradients are treated as independent uncertain variables that vary over a defined range, 
and the range is defined to be representative of an entire wind shear profile during descent. In this 
way, robustness over an entire wind shear profile can be explicitly incorporated into the design. It 
should be noted, however, that this approach should not "tune" the design to a particular wind 
profile since information about which gradient combinations will occur and in what sequence is not 
provided. Rather, this approach should provide robustness to all gradient combinations within a 
particular range of values, which are determined from a realistic but severe wind shear profile 
(such as the Kennedy microburst wind shear profile used in this study). Allowing the gradients to 
vary independently over all values within the designated range results in a rectangular uncertain 
parameter space. Intuitively it appears that some subspace (such as a diamond) within this 
parameter space might be a more accurate representation of the physical microburst wind shear 
phenomenon. However, further study of microburst wind shear data is required in order to 
substantiate that this is true in general. In addition, it is unclear whether some preferred subspace 
exists for realistically representing microbursts as well as other wind shear phenomena, such as 
terrain-induced and frontal wind shears, and what that subspace might be. Therefore, more 
extensive study into the structure of these various types of wind shear should be performed in 
order to define the optimal uncertain parameter subspace for wind shear, and to thereby avoid 
introducing unnecessary conservatism into the control system design. 

Actuator Model 

First-order approximations are used to represent the elevator dynamics, throttle dynamics, 
and engine dynamics. The state model associated with the actuator dynamics is given by the 
following equations. 



x A = A a x a + B a u c 

(11) 


u = C A x A 

(12) 

where: 

X A = [ 5th 5r r 

(13) 


u c = [ 8b c 5th c | T 

(14) 


and the matrices Aa, Ba, and Ca are defined in the Appendix. The control command inputs 
contained in the Uc vector are elevator position command (8^) and throttle position command 
(8thc)- These commands will be computed by the control law, K(s). 

An uncertainty model for this problem formulation is developed in the next subsection. 


2.2 Uncertainty Model Formulation 

Uncertainty can be introduced into a general system model as shown above in Figure 1. 

As stated earlier, A represents a block diagonal matrix of system uncertainties, and the additional 
Ry jfy and D matrices in G(s) represent the interconnection of the uncertainties into the system (see 
Belcastrcf, Chang, & Fischl (1991), ref 4). This uncertainty model will be formulated for the wind 
shear problem in this subsection. In this study, only uncertainty associated with the wind shear 
model was considered. In order to obtain the uncertainty model, the uncertain wind shear 
gradients shown in equation (10) must be separated from the nominal plant system. Details of the 
uncertainty model development are therefore given next. 
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2.2.1 Wind Shear Model Uncertainty 

As shown in equations (8-10), the A and B w matrices are linear (affine) matrix functions of 
the horizontal and vertical wind shear gradients U z and W z , respectively. These gradients were 
considered to be real uncertain variables which vary linearly with additive uncertainty, as follows: 

U z = U Z0 + 5u z (15) 

W z = Wzo + 5 Wz (16) 

where the nominal wind shear gradients, Uzo and W^, were set to zero (since wind shear 
penetration is not typical of a nominal landing), and the uncertain wind shear parameters, 8u z and 
5w z . were allowed to vary over a range defined by the wind shear gradients associated with the 
1975 Kennedy landing accident. In order to avoid existence problems in the control system design 
solution (as a result of setting Uzo W^ to zero), the integrator poles associated with the wind 
shear model are placed slightly into the left half-plane (e.g., -.001 as noted in the Appendix) 
during the design procedure so that stabilizability requirements for solution are met. The range of 
gradient values used for the uncertain wind shear parameters, 8u z and 8\v z , was determined using 
a simulation package developed by Stanford Research Institute (SRI) International, which 
reconstructs the wind profile that occurred during the 1975 Kennedy accident (see Dieudonne 
(1979), ref. [6]). For this wind shear profile, Uz and Wz were determined in Belcastro & Ostroff 
(1985), ref. [2], to vary as follows: 

(Increasing Headwind) -.25 < U z < .39 (Increasing Tailwind) 

(Increasing Updraft) -.47 <, W z < .33 (Increasing Downdraft) 

Therefore, for the purposes of this study, the wind shear uncertain parameters were allowed to 
vary as follows: 

I 8 Uz I < 0.4 , I 8 Wz I < 0.5 

where these ranges were selected to cover the wind gradient variation for the Kennedy wind profile 
described above. Nonsymmetric uncertainty bands, which are more representative of the actual 
ranges given above, will be included in future work. As described in Belcastro, Chang & Fischl 
(1991), ref. [3], Morton & McAfoos (1985), ref. |8|, and Manning & Banda (1989), ref. [7], the 
A and B w matrices can be expanded as follows: 

A = A 0 + Aa , B w = B W0 + B w a (17) 

and: M = [ Aa B wA j = M] 8u z + M2 8\v z (18) 

For this problem formulation, both Mi and M 2 are rank 2. Thus, as discussed in Morton & 
McAfoos (1985), ref. [8], and Manning & Banda (1989), ref. [7], the uncertain parameters 8u z 
and 8w z . will each have to be repeated in the uncertainty block, A, of the generalized uncertainty 
model (Figure 2). As discussed in Belcastro, Chang, & Fischl (1991), ref. [3], the uncertainty 
interconnection matrices B xp , D qw , C qx , D yp , and D qu matrices can be determined by decomposing 
the Mj matrices. However, since the uncertain parameters occur only in A and B w , the matrices 
Bxp* Dqw* and Cq X are nonzero and D yp , and D qu are zero. The Dm, matrix is zero because there 
are no crossterms of uncertain parameters (e.g., 8u z 5w z ) in any 01 the matrix elements. The B xp , 
D qw , Cn X , Dyp, D qu , and D qp matrices associated with this problem are given in the Appendix. 

In the generalized uncertainty model, it is desirable to scale the uncertain parameters so that: 

I8' Uz i< 1 , 18'wz'^l (19) 

in order to utilize the standard robust stability and performance |i-test given in Doyle, Wall & Stein 
(1982). Performing this normalization facilitates the analysis and design process in that the p 
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robustness test (without loss of generality) always involves comparing p. to a unity bound, as 
opposed to recalculating the bound for every problem - or even recalculating the bound within the 
same problem whenever there is a change in numerical values (such as the uncertainty ranges 
discussed here). This normalization was therefore accomplished by introducing the following 
scaling matrix, Asc: 


A sc = diag [0.4I 2 0.5 I 2 ] t (4x4) 

The A matrix associated with the uncertainty model is therefore given by: 


(20) 


A = diag [ 5 uj 2 SwJ 2 ] T (4x4) 


( 21 ) 


The manner in which A^ enters into the model will be shown in the summary of the general system 
model. The nominal plant model, P(s), consists of the matrices A 0 , B 0 , C 0 , D 0 , and D Wo , which 
are also given in the Appendix. 

2.2.2 Wind Turbulence Model Uncertainty 


For this problem formulation, no uncertainty was modeled for wind turbulence. However, 
the Dryden turbulence model (described in Dieudonne (1979), ref. [6], and included in the SRI 
wind simulation package) was used to provide turbulence disturbance weighting. This model will 
be presented later in Section 2.3.4. As will be discussed there, uncertain parameters associated 
with this model could be defined and incorporated into the uncertainty model. However, how to 
do this will be the subject of future work. 

2.3 Generalized System Model 


The next step in formulating the problem is to define the control and performance 
objectives, and obtain the associated weighting functions. This is accomplished in the following 
subsections. 


2.3.1 Control Objective 

As mentioned previously, the overall objective of this study is to track the glideslope, 
despite severe wind shear and turbulence conditions, during the descent to landing of the aircraft. 
The method used in this study to achieve this control objective was to minimize the error in h, V a , 
and 0 without exceeding actuator position and rate limits. Minimizing the error in h, which 
represents altitude along the glideslope, provides a means of minimizing deviations from the 
glideslope during descent, and minimizing errors in V a and 0 provides a means of minimizing the 
effects of wind shear on the aircraft dynamics. The controlled output vector, z, is therefore given 
by: 


z — C Zx x + C ZXA x A + D Z y u c + D/c r c (22) 

where: 

z = | eh e v , £e : 5 e 5,1, 6 C 8 U i 1 T = l Zc T Zc T ] T (23) 

r e = I h c 0 C 1 T (24) 

The first three elements in the controlled vector, z, represent errors associated with altitude, 
airspeed, and pitch attitude, respectively, and the last four elements represent elevator and throttle 
positions and rates, respectively. The matrices associated with these controlled variables are given 
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in the Appendix. The command vector, r c , represents the commanded altitude, airspeed, and pitch 
atti tu de , respectively. For this study, the control effort will focus on tracking a 3-degree 
g l id e, slo p r However, the problem was formulated so that glideslope capture and/or the flare 
maneuver to transition from the glideslope to touchdown could easily be implemented by 
incorporating appropriate command models for these variables associated with these maneuvers. 
For example, a command model for flare might include an exponential flight path command, he, 
with a linearly decreasing airspeed command, , and a linearly increasing pitch attitude 
command, 0 C . The equations associated with these commands would then be incorporated into the 
plant state model. 

2.3.2 Performance Objectives and Weighting Functions 

The performance objectives were formulated relative to the controlled variable vector, z, 
containing tracking errors and control command variables. The approach used in this study is 
similar to that described in Balas, Doyle, et al. (1991), ref. [1], and Bibel & Stalford (1991), ref. 
[4]. The performance requirements associated with glideslope tracking as well as the associated 
weighting functions are discussed first, followed by the control position and rate limitations of the 
control actuators along with their associated weights. 

The performance requirements associated with glideslope tracking error were specified 
relative to the following time response characteristics: 

Time Constant: 1.0 sec. 

Steady-State Error: 0. 1 % 

Step Response Overshoot: 10% 

The 1 second time constant corresponds to the natural frequency of the short period mode, and the 
other values were based on engineering judgement. It should be noted, however, that the design 
process is iterative, and that the above values are suggested as a reasonable starting point for doing 
a design for this vehicle. These requirements were then used to define a performance weighting 
function, Wp, of the following form: 


Wp(s) = 


Kp (s+coi) 
(s+0>2) 


The time constant requirement was used to establish the weighting function crossover frequency, 
tOc, i.e.: 


t = 1 .0 sec => CDc = 1 .0 rad/sec 

The steady-state error requirement was used to establish the dc gain of the step response signal 
(using the final value theorem) as follows: 


0.1% Steady-State Error => I W P (s) I l s _*o = 1/.001 = 1000 

where the dc gain is inverted in order to normalize the error response, thereby facilitating the 
analysis and design process (i.e., p can be related to unity in the standard p robustness test). 
Similarly, the percent overshoot requirement is used to specify transient response performance to a 
step response (using the initial value theorem): 

10% Overshoot => I W P (s) I l s -»oo = .1 

This approach is taken from the results presented in Bibel & Stalford (1991), ref. (4]. Using these 
specifications, W P (s) was determined as: 


8 



(25) 


W P (s) 


0.1 (s+10) 
(s+,001) 


and a frequency response plot is shown in Figure 3: 



0) (rad/sec) 


Figure 3. Frequency Response Log Magnitude Plot of the 
Performance Weighting Function, W p (s) 


This performance weighting function is related to the desired sensitivity function and is applied to 
the three error signals in z. However, if performance requirements for individual error signals 
vary, the above approach can be used to determine individual weighting functions for each. The 
state equations for the weighted error variables can then be written as follows: 


^w p — A Wp x Wp + B Wp Zq 
Zc = C Wp x Wp + D Wp Zc 


(26) 

(27) 


where Ze is defined using equations (22) and (23) and the matrix definitions given in the Appendix 
as follows: 


Zc = C/, X - 1 3 r c (28) 

and the state vector x Wp consists of the states associated with the three error signals in z being 
passed though the weighting function W p (s). The matrices A Wp , B Wp , C Wp , and D Wp in equations 
(26) and (27) are defined in the Appendix. 

The control position and rate signals can be weighted with constant scalars based on the 
actuator position and rate limits, as follows: 


Elevator Position: 
Throttle Position: 


6 c ^±10deg Trim = 3.2 deg 
10 < 5ih ^ 60 deg Trim = 17.6 deg 


Elevator Rate Li m i t: 10 deg/sec 

Throttle Rate Limit: 10 deg/sec 


I ^emax I ~ 2 deg 
I ^thmax ^ = 40 deg 
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Then the weightings on the control positions and rates were chosen to be: 

Wfie = 1/7, Wgih = 1/40, W&;(ratc) = 1/10, Wgihfraie) = 1/10 (29) 

where the limit values are inverted to normalize the controlled variable outputs - again, to facilitate 

the analysis and design process. 

The overall weighting matrix, W z (s), is therefore given by: 

W z (s) = diag [ W c W c 1 (30) 

where: W e = diag[Wp(s) Wp(s) Wp(s) ] (31) 

Wc = diag [ Wge W gc(ratc) ^5th(ratc) I (32) 

Obviously, this weighting matrix is varied during the design process to improve performance. It 
may even become necessary to vary the performance weighting function on the error signals 
independently. This weighting strategy is presented as a starting point for the design process. 

2.3.3 Disturbance Weighting 

The exogenous input weighting matrix, W w (s), is used to define the frequency spectra 
assoripM with the input disturbance vector, d w , measurement noise vector, r\, and command 
signal vector, re. For this problem formulation, W w (s) was defined as follows: 

W w (s) = diag | W g (s) I I ] (33) 

where Wg(s) is the weighting matrix associated with the wind gust disturbance vector, d w , defined 
in equation (5). The Etyden turbulence model given in Dieudonne (1979), ref. [6], was used to 
provide lowpass turbulence weighting functions for the longitudinal and vertical wind gust 
components, u wg and w wg , and this model is repeated below in Figure 4 for convenience: 


Longitudinal : 


Vertical: 


F u (s) 

U Wg 

- h 



Ww. 

F w (s) 

ft 

. m 




Figure 4. Dryden Turbulence Model 


In the above turbulence model n u and n w are white, zero mean, unit variance noise signals, V a 
represents trim airspeed, and the terms L u , L w , a u , and o w are defined for the Kennedy wind 
profile as follows in Table 1 : 
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Table 1. Turbulence Specifications for the Kennedy Wind Profile 

The output signals, u W g and w Wg , represent the longitudinal and vertical wind gust components 
comprising the d w vector, and these components are given in the body axes with units of ft/sec. 
The wind disturbance weighting matrix, W g (s), is therefore determined by the Dryden turbulence 
model as follows: 

W g (s) = diag [ F u (s) F w (s) ] (34) 

and the state model associated with this turbulence representation can be formulated as follows: 


x Wg - A w g x Wg + Dw g n Wg ( 35 ) 

y Wg = C Wg x Wg (35) 

where: Xw « = t Uw * Ww »‘ Ww * 2 ] T (37) 

n w g = [ n u n w ] T (38) 

y Wg = [ u Wg w Wg ] T ( 39) 


and the matrices A w „, B wg , and C wg are defined in the Appendix. The values from Table 1 that 
were used in the weighting functions were chosen to yield the highest bandwidth. These values 
correspond to the first row of the table. A frequency response plot of the longitudinal and vertical 
transfer functions at these values is shown below in Figure 5: 



to (rad/sec) 

Figure 5a. Frequency Response Plot of the Dryden Turbulence 
Model Using the Data Given in Row 1 of Table 1 
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co (rad/sec) 

Figure 5b. Frequency Response Plot of the Dryden Turbulence 
Model Using the Data Given in Row 1 of Table 1 


As mentioned previously, the Dryden model could also be used to develop an uncertainty 
model associated with wind turbulence, where the terms L u , L w , a u , and o w are considered to be 
real uncertain parameters which vary over the values defined in Table 1 for the Kennedy wind 
profile. As can be seen by the above transfer functions as well as the A wg , B W g, and C wg matrices 
given in the Appendix, these parameters enter into the model in a fairly complicated manner. Thus, 
modeling these parameters as uncertainties will be addressed in future work. 

As indicated in equation (33), no dynamic weighting was defined for the measurement 
noise vector, t|, or the command signal vector, r c . However, measurement noise weighting could 
be accounted for in the design process by defining an appropriate high-frequency noise spectrum 
for each sensor measurement signal and incorporating the associated states into the generalized 
model as demonstrated above for the turbulence model. Similarly, command signal prefiltering 
could be accounted for in the design model by defining appropriate weighting functions for the 
command signals in r c . 

2.3.4 Summary of the Generalized System Design Model 

The generalized system model as depicted in Figure 1 is now summarized. The vector 
states, inputs, and outputs for this problem are given as follows: 


X G = | X X A x Wg x w „ | T 

(40) 

P = 1 Pu Zl Pu Z2 Pw Zl Pw Z2 J T 

(41) 

w = [ dw T\ r c l T 

(42) 

q = [qu Zl qu Z2 qw Zl qw Z2 l T 

(43) 


The vectors x, xa, and x Wg are defined in equations (3), (13), and (37), respectively, and the x Wp 
vector represents the three states associated with the performance weighting functions for the three 
error signals in the control vector, z. The vectors u, y and z of Figure 1 are defined in equations 
(4), (6), and (23), respectively. The elements of d w , T), and r c are defined in (5), (7), and (24), 
respectively. The generalized plant model, G(s), contains the nominal plant model as well as all 
weighting and scaling matrices depicted in Figure 2 and used in the design and analysis. This 
model is summarized for the configuration of Figures 1 and 2 as follows: 
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The matrices A* B 0 , B Wq , C g , D 0 , and D Wo are associated with the nominal plant model, and the 
. D p r\ n anH D « are asscxiititcdi with the interconnection ot tne 

- associated whh ,he con.ro, >ed vanabie 

model, and the matrices A Wp , Bwp. Cwp, and D Wf are associated I with the performance 
weighting matrix, W e , associated with the performance weighting function Wp(s). All ot these 
matrices are defined in the Appendix. The uncertain parameter scaling matrix, A*, 18 “ 

equation (20), and the control weighting matrix, W c , is defined by equations (32) and (29). The A 
Sx for this problem is defined in equation (21 ). The only other component m the generalized 
model of Figure 1 is the controller, K(s), which could now be designed using p-synthesis. The 
controller design for this problem formulation will be the subject of future work. 
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3. Conclusions 


This paper has presented a formulation of a glideslope tracking problem through severe 
wind shear and gusts for application to robust control system design using p-synthesis. The 
uncertainty in this problem is associated with the wind shear model and appears as real parameter 
variations of the horizontal and vertical wind shear gradients. The range of parameter variation 
was based on the gradient extremes determined from a reconstruction of the wind profile that 
occurred in the 1975 landing accident at Kennedy International Airport. Formulating this problem 
for p-synthesis should provide a means of incorporating robustness to wind shear explicitly into 
the design process without "tuning" the design to a particular wind gradient or profile. Wind gusts 
were also included in the model as disturbance inputs to the plant. The Dryden turbulence model 
was used to provide weighting information. However, this model could also be used to formulate 
additional wind model uncertainty by treating the turbulence parameters L u> L w , o u , and o w as real 
uncertain parameters. Future work in the area of wind shear and gust uncertainty modeling will 
address the inclusion of these turbulence terms as part of the uncertainty model, as well as the 
determination of an appropriate wind shear parameter subspace. Nonsymmetric uncertainty bands 
associated with the uncertain parameters will also be addressed. Performance weighting functions 
were also derived for inclusion in the problem formulation developed in this paper. This problem 
formulation can be used to obtain a p-synthesis control system design, although modifications to 
the design parameters (e.g., the weighting functions) may be required to obtain a final design. The 
results of the control system design as well as a closed-loop evaluation using a full nonlinear 
simulation of the aircraft will be presented in a subsequent paper. 
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Appendix: Model Definition 



where: 

‘ -.41790E-01 .98126E-01 -.49565E+01 -.32155E+02 

A _ -.29344E+00 -.76071E+00 .21178E+03 .96749E+00 

“ " -.30862E-03 -.49839E-02 -.52858E+00 -.37077E-03 

. .00000E+00 .00000E-KX) .1(XXX)E+01 .00000E+00 - 

' .41790E-01 -.98126E-01 

A _ .29344E+00 .76071E+00 A _ [ a w ii a wJ2 a w n a w i 4 ' 

* w .30862E-03 .49839E-02 ’ wa l a w2 i a w22 a w23 a W 24 . 

. .(XXXX)E+0() ,(X)(XX)E+(X) - 

a w n = .028944 U z + .0(X)83845 W z , a w2 i = -.00083845 U z + .028944 W z 
a w i 2 = .99916 U z + .028944 W z , a w 22 = -.028944 U z + .99916 W z 

a w i 3 — 0 > a w 23 = 0 

a w i 4 = -210.89 U z - 17.188 W z , a w2 4 = 17.188 U z - 210.89 W z 

U z = Horizontal Wind Shear Gradient (Uncertain Parameter) 

W z = Vertical Wind Shear Gradient (Uncertain Parameter) 


A 


WW 


a w i5 0 
0 a w2 6 . 


aw 15 = 0 (aw 15 = — .(X) 1 for Design Model) 
aw 16 = 0 (aw 16 = -.001 for Design Model) 


A h = [-.028956 -.99958 0 211.3] 


.40382E-02 .40218E-03 

n - [ B «l r - -.17254E+00 -.20149E-06 

[OJ 1 J ’ a -.21025E-01 .62988E-05 

.0(XX)0E+(X) .OOOOOE+OO 
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c = 


.9997 .0234 0 0 -.9997 -.0234 0 " 

0 0 10 0 0 0 

0 0 0 1 0 0 0 

0 0 0 0 0 0 1 

-.028956 -.99958 0 211.3 0 0 0 

D = [0] (5x2) , D w = [0] (5x2) 


Plant Trim Values: 


Variable 

Trim Value 

Units 

u 

211.53575 

ft/sec 

w 

4.9508090 

ft/sec 

q 

0 

rad/sec 

6 

-.028960023 

rad 

U w 

0 

ft/sec 

w w 

0 

ft/sec 

h 

200 

ft 

6c 

3.178 

deg 

6th 

17.576 

deg 

8t 

8,112.717 

lbs 

V a 

211.5937 

ft/sec 

U z 

0 

ft/sec/ft 

W z 

0 

ft/sec/ft 


Actuator State Model 



-W5. 

0 

0 



0 1 

A a = 

0 

0 

-<qs th 

^CngflSen, 

0 

. B a = 

0 

0 

1 

3° 


where: 


( 05 e = 16 rad/sec, G)s t h = 10 rad/sec, wgcng = 0.5 rad/sec, Kc ng = 596 lbs/deg 


C A 


1 0 0 
0 0 1 


Uncertainty Model 


Nominal Plant Model: 

Aq = A (A wa — 0) , B 0 = B , Cq = C , Do = D, Dw 0 = Dw 
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Uncertainty Matrices: 


B X p - 


0 0 0 0 

6 6 6 6 
1010 
0 10 1 
oooo 


(7x4) , D qw = 


-.028944 -.99916 

.00083845 .028944 
-.00083845 -.028944 
-.028944 -.99916 


Cq X - 


.028944 .99916 0 -210.89 0 ••• 0 

-.00083845 -.028944 0 17.188 0 ••• 0 
.00083845 .028944 0 -17.188 0 0 

.028944 .99916 0 -210.89 0 0 


(4x7) 


D yp = f0| (5x4) , Dq U = |0] (4x2) , D qp = [0] (4x4) 


Controlled Variable Model 

c, - [ c„ C,„ ] = 


C-U 

0 


0 

C*. 





10 0' 


0 0000 01 


0 1 0 

= 

.9997 .0234 0 0 -.9997 -.0234 0 

> C /tt - 

-0)6 e 0 0 


0 0 0 1 0 0 0- 


0 -^u, 0 - 


D*. = 


0 


D- 


Zuu 


(7x2) , D y<lll = 


D/.= 


0 


0 

0 

w 5. 

0 


(7x3) 


0 

0 

0 


Performance Weight ing Model 


Aw P = 


-.001 

0 

0 


0.1 

0 

0 

0 

-.001 

0 

, Bw P = 

0 

0.1 

0 

. 0 

0 

-. 001 . 


. 0 

0 

0.1 . 


Cw P = 


■ 9.999 

0 

0 


' 0.1 

0 

0 

0 

9.999 

0 

, Dwp = 

0 

0.1 

0 

. 0 

0 

9.999 - 


- 0 

0 

0.1 . 
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Wind Turbulence.. Model 



Va 

U 

0 

0 



0 

0 

Va 

Lw 




0 


-V, 







10 0 
0 0 1 . 
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